Dynamics of the two-dimensional 5=1/2 dimer system (CsHe^F^CuCU 
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Inelastic neutron scattering was used to study a quantum 5 = 1/2 antiferromagnetic Heisenberg 
system — Bis(2-amino-5-nuoropyridinium) Tetrachlorocuprate(II). The magnetic excitation spec- 
trum was shown to be dominated by long-lived excitations with an energy gap of A= 1.07(3) 
meV. The measured dispersion relation is consistent with a simple two-dimensional square lattice 
of weakly- coupled spin dimers. Comparing the data to a random phase approximation treatment of 
this model gives the intra-dimer and inter-dimer exchange constants J = 1.45(2) meV and J / =0.31(3) 
meV, respectively. 



PACS numbers: 75.10.Jm, 75.50.Ee 

Quantum spin liquids, where strong zero point fluctua- 
tion destroy magnetic long-range order even at T=0, have 
been a subject of intense studies for over two decades. Of 
particular interest are two-dimensional systems of this 
type. In addition to their relevance to layered supercon- 
ducting cuprates^ they are of interest in the context 
of field-induced quantum phase transitions^ where, to 
date, most experimental work was performed on either 
spin chains^ and ladders^ or 3D systems i 10 ! 11 Several 
2D spin liquids including BaCuSi 2 6 r^ SrCu 2 (B0 3 )2r 3 - 
and PHCCi^ have been studied in great detail, but show 
rather complex behavior due to geometrically frustrated 
interactions. 

Recently, a new quantum antiferromagnet Bis(2- 
amino-5-fluoropyridinium) Tetrachlorocuprate(II) 
(5FAP-CuCi4, for short) has been discovered and 
characterized by bulk measurements^ The ground 
state is a spin liquid. The energy gap A to the lowest 
excited (triplet) state was determined to be 0.97(4) 
meV from the high-field magnetization data. 5FAP- 
CuCU ((CsHe^F^CuCU) crystallizes in a monoclinic 
space group P2\/c with lattice constants a=6. 926(7) 
A, 6=21.73(2) A, c=10.911(10) A, and /3=100.19(2)°. 
There are four inequivalent Cu atoms per unit cell as 
shown in Fig. Ufa). In view of crystal structure, 5FAP- 
CuCU consists of CuCl^ - anion layers separated by 
organic cations. Those CuCl^ - anions form zigzag layers 
which are stacked along the the crystallographic a axis 
as shown in Fig. [2(b). Fig. [TJc) shows the arrangement 
of CuCl^ - anions and organic cations viewing from the 
a axis. The yellow and purple lines indicate the most 
plausible super-exchange interactions that could be 
carried across Cu-Cl-Cl-Cu bonds. The former involves 
the shortest Cl-Cl contact with a separation of 3.66 
A, while for the latter the relevant Cl-Cl distance is 
4.07 A. One can expect inter-layer coupling along the 



a axis to be exponentially weak due to a much larger 
Cl-Cl separation 5.07 A. Even with this in mind, it 
is not obvious which in-lane interactions are the most 
relevant. To better understand the spin Hamiltonian 
for 5FAP-CuCl4, in this paper, we report inelastic 
neutron scattering (INS) measurements of magnetic 
excitations and their dispersion relation. We show 
that 5FAP-CuCU can be adequately described by a 
minimal model that is equivalent to a square lattice of 
weakly interacting antiferromagnetic S — 1/2 dimers 
with no geometric frustration. 

Small size fully deuterated single crystalline samples 
of 5FAP-CuCl4 were prepared by the procedure described 
in Ref. [l5|. About twenty single crystals (~0.3 g to- 
tal) were co-aligned with a mosaic about 4° and used 
for INS measurements on the cold-neutron triple-axis 
spectrometer IN14 at Institut Laue-Langevin. Another 
assembly consisting of 20 single crystals (~0.5 g) with 
a total 5° mosaic was used for the measurements on 
TASP — spectrometer at SINQ, Paul Scherrer Institute. 
Horizontal beam divergences were given by 58 Ni guide- 
60'-open-open at IN14 and by lO'/A-open-open-open at 
TASP. Both experiments were performed with the scat- 
tered wave vector kf fixed at 1.5 A -1 and Be filter placed 
after the sample was used to remove higher order beam 
contamination. The sample was oriented horizontally in 
a (0,fc,Z) reciprocal lattice plane and cooled in a standard 
He-4 cryostat. 

The data were collected in a series of constant- q scans 
at T = 1.5 K along the high symmetry directions in the 
(0,fc,Z) reciprocal plane. The inset of Fig. [2fa) shows the 
raw data of a typical scan performed at g=(0,0,-1.5), 
where a sharp resolution-limited peak due to a magnetic 
excitation is clearly visible. The background (shown as 
a solid line) was determined by fitting the data over the 
range where no magnetic excitations are apparent, to a 
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FIG. 1: (color online) (a) Outline of a unit cell in 5FAP- 
CuCU. (b) View from the c axis, showing the zigzag layer 
structure. Only Cu and CI atoms are shown, the organic 
cations have been removed for clarity, (c) View from the a 
axis, showing the arrangement of CuCl^ - anions and organic 
cations. Different color lines stand for the different magnetic 
interactions. Color coding is as follows: Cu, CI, C, H, F, and 
N. 



term linear in energy, plus a resolution-limited Gaussian 
peak at zero energy transfer. The latter accounted for in- 
coherent elastic scattering. Background-subtracted data 
at several wave vectors are shown in Fig. 02 At certain 
wave vectors the observed signal has a distinct double- 
peak structure as shown in Figs.[3fc) and (d). Solid lines 
in Fig. [3] are fits to Gaussian profiles of fixed widths, 
determined by instrumental resolution, as calculated us- 
ing the Cooper-Nathans approximation^ Excitation en- 
ergies (dispersion relation) obtained from such fits are 
summarized in Fig. [2j The spectrum appears to consist 
of two equivalent branches, an "optic" and an "acoustic" 
one, although the former is not visible at the majority of 
the wave vectors studied. 

The most comprehensive method of determining the 
relevant exchange parameters in a gapped system with 
weakly dispersive excitations is based on the application 
of the first-moment sum rule for the magnetic dynamic 
structure factor J£ This procedure involves quantitative 
measurements of inelastic intensities as a function of wave 
vector and was successfully applied to the rather complex 
interaction geometry in CuHpCl. 10 Due to technical com- 
plications, associated with the sample fracturing upon 
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FIG. 2: Dispersion of magnetic excitations measured in 5FAP- 
CuCU at T = 1.5 K as a function of k and /. The solid 
(dashed) line is a fit to the acoustic (optical) magnon dis- 
persion mode as discussed in the text. Inset: Constant- 
q = (0,0, —1.5) scan as a function of the transferred energy 
fvjj at T = 1.5 K. The solid line is a Gaussian and a linear 
term fit to account for the background contribution. 



cooling, and the data being collected in several disjoint 
experiments, in the present study we adopted a simpler 
approach. Only the measured excitation energies were 
analyzed. A parameterized model was refined to repro- 
duce the experimental dispersion relations. The proposed 
minimal model is visualized in Fig. |4j which shows the 
exchange constants a single layer of magnetic Cu 2+ ions. 
The dimers are formed by the stronger antiferromagnetic 
bonds J, which ensures a spin-singlet ground state and 
an energy gap. The dispersion of singlet- triplet dimer ex- 
citations is due to the weaker inter-dimer exchange con- 
stant J' . Topologically, this spin network can be viewed 
as a square lattice of S = 1/2 dimers. Nevertheless, its 
actual realization in the crystal has two dimers per crys- 
tallographic unit cell. This leads to a folding of the Bril- 
louin zone and the appearance of two excitation branches. 
As shown below, this minimal model seems to be fully 
consistent with experimental observations. 

We can treat the weak inter-dimer coupling J' at 
the Random Phase Approximation (RPA) level. 19 This 
method has been applied to numerous spin systems, 
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FIG. 3: Representative background-subtracted INS data 
for 5FAP-CuCl 4 at T = 1.5 K. The solid lines are fits Gaus- 
sian profiles of instrumental width. The shaded area in (a) is 
excluded due to a contamination by "accidental Bragg" spu- 
rious scattering. 



cell as the basic spin cluster (tetramer). These tetramers 
will then form a Bravais lattice. With four sites in each 
tetramer, their bare susceptibilities are written as a 4 x 4 
matrix, the indexes enumerating the four spins in each 
unit: 
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is the difference in popula- 



te R ( T ) = We -JHk B T) 

tion of the excited and ground states of the dimers. The 
weak interaction J' couples adjacent tetramers, but also 
dimers within each tetramer. Its effect is taken into ac- 
count within the RPA. The corresponding RPA equation 
is written as: 



x RPA M = xoH[i - xoMJ'(q)]- 1 . 



(2) 



Here J'(q) is the Fourier transform of exchange interac- 
tions due to J' . It is also written as a 4 x 4 matrix: 



J'(q) = J'l 
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FIG. 4: (color online) Minimal model for magnetic inter- 
actions in 5FAP-CuCl4, shown in projection onto the (6, c) 
plane. Dimer singlets are formed by antiferromagnetic ex- 
change J. Dispersion of singlet-triplet excitations is due to 
inter-dimer interaction J' . 
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including Cs 3 Cr 2 Br 9f 20 BaCuSi 2 6 , 21 BaCuSi 2 7 
Ba 3 Mn 2 8 r^ and Ba 3 Cr 2 8 ^ In application to 5FAP- 
CuCU , we can treat two decoupled dimers in a unit 



The dispersion relation is then readily obtained from 
the location of poles in x RPA { Q ii U0 )\ 

(hu^) 2 = J 2 =F JJ'R(T) [costt(/c + I) + costt(/c - I)] . 

(4) 

The solid and dashed lines in Fig. [2] are fits of this 
dispersion relation to the experimental data. A good 
agreement is obtained with J = 1.45(2) meV and J' = 
0.31(3) meV. 

In summary, 5FAP-CuCl4 can be described as a simple 
2D network of weakly interacting spin dimers, which is 
topologically equivalent to a dimer square lattice. The 
magnon dispersion bandwidth is 0.70(4) meV, smaller 
than gap energy A = 1.07(3) meV, but still substantial. 
Future work will focus on the high-field study of quantum 
phase transition and critical phenomena in 5FAP-CUCI4. 
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